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Abstract 

We discuss effects of particle interaction on Bose condensation in inhomogeneous traps with and 
witliout optical lattice. Interaction pushes normal particles away from the condensate droplet, 
which is located in the center of the trap, towards the periphery of the trap where the trapping 
potential is large. In the end, the remaining normal particles are squeezed to a quasi-2D shell 
around the condensate droplet thus changing the effective dimensionality of the system. In the 
absence of the optical lattice the index in the temperature dependence of the condensate density 
at the later stages of the process is close to 2 with a weak dependence on the number of trapped 
particles. In the presence of the lattice inside the trap this index acquires a strong dependence 
on the number of particles inside the trap and gradually falls from a 3D to a 2D value with an 
increase in the number of particles. This change in index is explained by the lattice-driven spread 
of the condensate droplet and the localization of the narrow band particles by the trap potential. 
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le interplay between the repulsive interac- 
^ since these two factors have opposite 



The spectacular experimental discovery of Bose condensation (BEC) made the study of al- 
kali gases in traps the focal point in atomic, low temperature, and condensed matter physics. 

For the first time, it became possible to observe some of the phenomena that have been dis- 

P 

cussed earlier only within theoretical models (see review |1]). The phenomena in ultracold 
alkali gases are incredibly rich and combine features inherent to diverse condensed matter 
and low temperature systems (Refs. 2|] and references therein) from "classical" superfiuid 
or superconducting systems to spin-polarized quantum gases ^ to Mott transition in the 
optical lattice 

The trap potential is always inhomogeneous. T : 
tion and the trapping potential complicates BEC 
effects on condensation: while the trap tends to concentrate the condensate in a narrow 
region of space around the particle ground state in the trap, the repulsion is responsible for 
the widening of this condensate droplet. The analytical description of the combined effects 
tends to be rather elusive and our previous experience with condensation in homogeneous 
systems is not very helpful. The problem becomes even more complex in the presence of 
the optical lattice inside the trap which adds two different localization processes - Mott 
transition and localization of narrow band particles by an inhomogeneous potential. 

Below we investigate a situation in which it is possible to get an accurate picture of the 
condensation in trapped interacting gases. The main attention is paid to the index in the 
temperature dependence of the condensate fraction and to the size of the condensate droplet. 
It turns out that this index is not universal even for a low density gas. What is more, the 
effective dimensionality of the problem changes with condensation and the later stages of 
BEC are different from initial. 

We start from BEC in trapped gases without the optical lattice, and add the complications 
associated with the optical lattice later on. We assume that the density is sufficiently low 
to neglect the interaction before the onset of condensation even in the center of the trap. 
This means that Tc is unaffected by the interaction. The interaction is brought into play 
only with the onset of condensation since the particles condensate in the center of the trap 
making the density in the center large. This makes the interaction, which is proportional 
to the particle density, large only in and around the condensate droplet. This also means 
that the normal particles are pushed out by the dense condensate towards the periphery 
of the trap where the interaction is negligible. The further particles move away from the 
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center the higher is the gradient of the trapping potential which is responsible for the force 
pushing the normal particles back towards the trap center. Thus, at the later stages of 
BEC, the majority of remaining normal particles are located in an almost two-dimensional 
shell around the condensate droplet and the dimensionality of the problem changes from the 
three-dimensional in the beginning of the condensation to quasi-2D later on. 

We consider a 3D harmonic trap with a single-particle ground state of frequency uj and 
spacial size cxo (axial asjmametry of real traps is irrelevant in our context). Without inter- 
action, BEC starts at Tc = [7] and the size of the condensate droplet is (Tq. 
Repulsion increases the size of the droplet with Nc (T) particles to a (T). Then the potential 
well for normal particles U (r) has a shell-type structure. 



t;(r,r) = ifiv 



where No = {^/^T/8) umaQ/has and we assume that the condensate density is Gaussian. 
The number of normal particles Nn (T) = N ~ Nc (T) is determined from the condition 
fi = 0. The size of the condensate droplet a (T) can be obtained from minimization of 
the condensate energy, including repulsion, similarly to Ref. js]. The interaction between 
the normal particles can often be excluded from Eq. ([T]). First, for less than 10^ particles 
in a trap, the density of the normal particles is negligible even in the trap center. For 
larger A^, the number of the normal particles on the later stages of the condensation is 
small. Finally, the density of the normal particles is suppressed even more by repulsion from 
the condensate droplet which spreads them through a large shell around the droplet Ana'^aQ 
instead of concentrating them near the center in the volume (47r/3) ctq. This gives at least 
an extra order of magnitude for which we can neglect the interaction of normal particles. 

No in Eq. ([T]) is the minimal number of particles in the condensate that is sufficient to 
create a strong repulsive core in the center of the trap. When N > Nc ^ No the normal 
particles are pushed away from the center by the repulsive core ([1]) into a potential valley 
surrounding the condensate droplet. For Rb in a trap with uj = 24 Hz, the values = 58.2 
A, (To = 2.2 X 10^^ m, and the critical number A^o that changes the topology of the normal 
cloud is A^o ~ 84. The center of the trap becomes inaccessible for normal particles when T 
is much smaller than the repulsion from the core. Using Tc instead of T and A^ instead of 
Nc, one gets a^No/a^ <^ A^^/s ^^^^^ ^^le critical value of Nc is around 10^. All this means 
that our results are applicable for A^ in the range 10^ -r- 10^. 



We are able to obtain an analytical description of the situation [cf. Ref. Gj). At the later 
stages of the condensation, the potential ([1]) forms a distinct valley away from the center of 
the trap as soon as N^. ^ Nq and equations for Nc (T) and a (T) reduce to 
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with 13 = huj/T, X = a/2 In (2x). The summation provides the temperature dependencies 
N,{T) and a(T). 

We found that the condensate fraction at the later stages of condensation can be given 

as 



NjN = 1 - (t/t;)" (3) 

with a relatively high accuracy. The important feature of Eq.(l3]) is that the temperature 
is normalized not by the critical temperature for the onset of condensation but by a 
different value T*. Since the squeezing of the normal particles towards the fringes of the 
trap accelerates with the number of particles in the condensate the normal shell narrows 
with increasing Nc, and, therefore, A^. As a result, the effective temperature T* should 
be higher than Tc and increase with increasing A^. Dependence of T*, or, more precisely, 
T* /hioN^^^ , on A^ is presented in Figure 1. For comparison, the critical temperature Tc for 
non-interacting particles in a 3D harmonic trap is Tc = 0.9hujN^^^ [t]. 

The striking change in behavior of T* (N) in Figure 1 occurs at A^ for which Tc ~ 
^Jtlj {NcCt'q / Noa^) . At higher densities the repulsion from the condensate droplet keeps 
the normal particles near the bottom of the potential valley around the droplet; at lower 
densities, the normal particles spread out and can even reach the center of the trap. An 
anomaly at the same threshold density is also observed in a (N), Figure 2, though the index 
a remains very close to the value 2 and is practically independent of A^, a = 2.02 ± 1%, in a 
wide range of A^ from 10'^ to 10^. This weak dependence a (N) is surprising for a nonlinear 
problem of this nature. The residual temperature dependence a (T) is within the same error 
bars. 
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FIG. 1: Density dependence of the reduced critical temperature T* /hujN^^'^ , Eq. 
interacting gas in a 3D harmonic trap this ratio should be 0.91. 
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FIG. 2: Density dependence of the index a, Eq.([3]) . For a non-interacting gas in a 3-D harmonic 
trap this index should be 3. 



These results confirm the evolution of the effective dimensionality from 3D, for which 
a = 3, to quasi-2D and the effective narrowing of the trap during condensation. 

The situation with an optical lattice (Refs. and references therein) with a period ao 
inside the trap is more complex. Here one deals with the Hubbard Hamiltonian, modified by 
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the trap potential, and can encounter the Mott transition jol which requires full occupancy 
of the lattice sites. The latter can occur with lowering of the temperature when particles 
gravitate towards the bottom (center) of the trap. With sufficiently strong on-site repulsion, 
the localization is practically inevitable for the condensate in the center of the trap though, 
of course, the Mott transition is sensitive to the trap profile jo], [l^ . The increased size of 
the condensate droplet in comparison to the system without the lattice changes the normal 
cloud surrounding the condensate for which it is possible to disregard the Mott transition. 

In the case of low initial density of particles na^ <^ 1 and strong on-site repulsion, the 
condensation starts at the same temperature as without the interaction. The condensate 
forms in the center of the trap and rapidly expands in size because of the strong on-site 
repulsion which tends to keep the density rical ^1. As a result, the size of the condensate 

1 /3 

droplet a ~ a^Nc becomes larger than o"max ~ (2 ^ 5) ctq for traps without the optical 
lattice. We will not dwell on potential freezing of the condensate resulting from the Mott 
transition and will concentrate on the condensation of the normal gas outside the condensate 
droplet. 

The main changes are associated with the band nature of the energy spectrum for particles 
in the optical lattice and a more complicated form of the wave functions. For the sake of 
comparison, in computations we use the same trap potential and particle scattering length. 
For the particle effective mass we use in computations the value 0, m* = 16m. 

The single-particle spectrum in the optical lattice e (p) has a band structure with a band- 
width A. The effect of the trap potential Utr (r) = ^hu {r/aof on the particles with narrow 
bands results in localization of particles with energy E in 2D shells e (p) + Utr {r) = E of the 
thickness £ (r) ~ {A/hu) cTq/t. [An exception is the center of the trap, where the gradient of 
the potential is small] . The particle wave function consists of three regions: rapid oscillations 
within this classically accessible shell and two attenuating tails beyond the classical turning 
points. The wave function for a particle with the energy E decays relatively slowly beyond 
the turning points, often as the Airy function of the type 3;+ (m*A/4/i^z/^)^^'' j 

where x is the distance from the "center" of the classically accessible shell for the particle 
with the energy E in the direction of the gradient, u = ^J2hjj [E — A/2)/cro, and m* is the 
particle effective mass at the turning point. The spatial distribution of particles should be 
calculated taking into account all three regions since for relatively shallow traps the contri- 
bution from the tails of the wave function can be large. Since such localization suppresses 
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FIG. 3: (Color online) Size of the condensate droplet a {T = 0) relative to the size of the trap 
(To, cr/cTQ, with (curve 1) and without (curve 2) the optical lattice as a function of the number 
of particles in the trap. The scattering lengths and effective masses are identical in both cases. 
Parameters of the lattice and the trap are given in the text. 

the accessibility range of narrow-band particles, the density in each point contains the con- 
tributions from the particles in a finite range of energies that are localized close to this point. 
For example, since only the particles with very low energies, E < A, can reach the center 
of the trap, the density in the center is suppressed in comparison with the trap without the 
optical lattice inside. 

As above, we start from the situation when the particle density above condensation is 
low and the (Hubbard) repulsion in the normal phase is negligible. The condition of low 
density allows us also to disregard the Mott transition in the normal phase Since the 
particles in the optical lattice are located mostly on the lattice sites of the size rather 
than spread uniformly, the repulsion is more effective than without the lattice. This means 
that the size of the condensate droplet a (T) should be larger than in the absence of the 
lattice. This is illustrated in Figure 3 which presents the ratio a (T = 0) /ctq for identical 
traps with (curve 1) and without (curve 2) the optical lattice. The scattering amplitude a^, 
which is responsible for repulsion, is the same in both cases. 

This seemingly innocuous lattice-driven increase in a leads to major effects and can 
eliminate a repulsive bump ([T]) in the center of the trap (at Nca^/NQa^ = 1) thus restoring 
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FIG. 4: Index a, Eq.Q, as a function of the number of the trapped particles in the presence of 
the optical lattice. Parameters of the lattice used in the computation are given in the text. 

the potential's parabolic structure in the central area. As a result, the change in index a, 
Eq. ([3]), is even more dramatic than the change in a, Figure 4. 

In Fig. 4, a starts from a 3D value at small density of particles which is understandable 
since there is no repulsive core in the center. With increasing number of particles the size of 
the condensate droplet grows leaving fewer normal particles in the central area and gradually 
reducing a to its quasi-2D value. What is not clear is why does a continue to decline with 
a further increase in A^; however, since our approach loses accuracy beyond = 10^, we do 
not present these data in the Figure. In general, the decrease in a (N) is accompanied by 
an increase in T* (N), which in the presence of the optical lattice grows much faster than 
A^^/^-dependence inherent to a free gas in a trap. 

In summary, we calculated the index for a temperature dependence of the condensate 
fraction for interacting gas inside harmonic trap. The results for traps without the optical 
lattice inside are quite clear: the repulsion from the condensate droplet pushes normal 
particles away from the center of the trap and concentrates them in a relatively thin shell 
around this droplet. Then the condensation becomes almost quasi-2D with the index a ~ 2. 
The presence of the optical lattice inside the trap changes the situation. The index a acquires 
a strong dependence on the number of particles inside the trap and gradually falls from a 
3D to a 2D value with an increase in the number of particles. This change in the index. 
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which is caused by the presence of the optical lattice, is explained by the wider spread of the 
condensate droplet and the localization of the narrow band particles by the trap potential. 
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